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ABSTRACT 


Generalized  voronoi  diagrams  of  sets  of  digital  curves  are 
a  helpful  tool  in  picture  analysis.  In  this  paper,  an  algori¬ 
thm  for  computing  labeled  Voronoi  diagrams  for  digital  straight 
line  segments  is  given.  Special  emphasis  was  given  to  the  use 
of  a  labeled  Euclidean  distance  transform.  This  transform  is 
the  key  feature  of  the  proposed  label  propagation  process.  The 
proposed  parallel  algorithm  for  computing  labeled  Voronoi  diagrams 
has  time  complexity  0(max{M,N})  for  input  pictures  of  size  Nxm 
using  a  mesh-connected  array  processor.  The  proposed  serial 
algorithm  for  computing  labeled  Voronoi  diagrams  has  time  com-  , 
plexity  0(MN).  ~  ^ 
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1.  Introduction 

At  some  stage  in  the  process  of  image  analysis,  a  given  image 
may  be  represented  as  a  set  of  points,  a  set  of  lines,  or  a  set 
of  curves.'  Clusters  of  these  sets  often  correspond  to  meeuiingful 
objects  in  the  scene.  In  this  paper,  we  introduce  a  labeled 
digital  Euclide^ul  distance  trauisform  for  constructing  labeled 
Voronoi  polygons  of  a  set  of  curve  segments. 

Consider  a  finite  set  P  of  points  in  the  pleme.  One  of  the 
basic  geometrical  structures  which  C2ui  be  used  to  characterize 
the  shape  of  such  a  set  P  is  a  partition  of  the  plane  into  a  set 
of  disjoint*  polygons  known  in  computational  geometry  as  the 
(nearest-point)  Voronoi  diagram  or  the  Dirichlet  tessellation 
[1] .  Here  each  polygon  consists  of  those  points  of  the  pleuie  that 
are  closer  to  some  given  point  of  P  tham  to  amy  other  point  of  P. 
The  Voronoi  diagram  depends  upon  the  particular  metric  that  was 
used. 

Voronoi  diagreun  structures  have  been  used  in  geography  [2]  as 
Thiessen  polygons,  in  physics  as  Wigner-Seitz  cells,  emd  also  as  a 
data  auialysis  tool  in  statistics  [3].  Recently  they  have  been 
used  as  an  alternative  way  of  extracting  skeletons  of  objects 
[4^5]  as  well  as  for  clustering  and  euialyzing  dot  patterns  in 
pattern  recognition  [6-8] . 

In  the  Euclidean  plane  with  the  Euclidean  metric  d,  the 
Voronoi  diagram  of  a  set  of  points  P*{pl,p2, . . . ,pn}  may  be 

*To  be  more  exact,  the  interior  sets  of  these  polygons  are 
disjoint,  but  different  polygons  may  share  an  edge.  Edges  of 
polygons  will  be  called  Voronoi  edges. 


represented  by  a  set  of  disjoint  convex  polygons  {G1,G2, . . . ,Gn}, 
where  a  point  x  of  the  real  pleuie  is  in  polygon  Gi  iff  there  is 
no  other  point  pj  in  P  which  is  closer  (according  to  metric  d)  to 
X  than  point  pi,  for  lsi£n.  Formally, 

Gi  ■  {x  :  d(x,pi)  *  Min{d(x,pj) ?  i<j<n}},  for  I5i<n. 

Although  Voronoi  diagrams  originally  were  designed  for  point 
sets  in  m-dimensional  Euclide2ui  space,  they  c2Ui  be  easily  ex¬ 
tended  to  sets  whose  elements  are  lines,  curves,  or  other  sim¬ 
ple  shapes.  Let  S  be  a  pleuiar  compact  point  set,  euid  let  p  be  a 
point  in  the  pl2me.  Then  d(p,S)  is  defined  to  be  the  minimal 
Euclidean  dist£Lnce  from  point  p  to  set  S,  i.e.,  d(p,S)»inf {d(p, 
q)  :  q€S}. 

The  Voronoi  diagram  for  a  set  {P1,P2 , . . . ,Pn}  of  lines  or 
curves  in  Euclideeui  space  can  be  defined  as  a  set  {G1,G2 , . . .  ,Gn} 
of  closed  point  sets  where 

Gi  =  {x  1  d(x,Pi)  =  inin{d(x,Pj)  :  isjsn}},  for  lsi<n. 

M^uly  algorithms  have  been  introduced  for  computing  voronoi 
diagrams  of  point  sets  in  two-dimensional  Euclidean  space  [9-131. 
The  diagrams  ceui  be  constructed  by  computing  one  polygon  at  a 
time  (simple  but  slow) ,  or  by  constructing  auid  merging  polygons 
using  the  divide-and-conquer  technique  which  results  in  algori¬ 
thms  with  time  conplexity  O(nlogn).  Recently,  O(nlogn)  algo¬ 
rithms  have  been  proposed  for  cosqputing  the  Voronoi  diagram  of 
a  set  of  line  segments  or  circles  [4,5]. 

Matsuyeuna  [14]  has  described  an  algorithm  for  extracting  the 
medial  axis  of  a  simple  polygon  from  the  Voronoi  diagreun  of 
disjoint  boundary  segments  of  the  object;  i.e.,  only  an  incomplete 


border  of  the  object  is  available  as  input.  (For  the  definition 
of  medial  axia  see  [15],e.g.;  roughly  speedcing^  this  is  the  set 
of  all  centers  of  nuocimal  disks  contained  in  the  given  object.) 

The  al^rithm  in  [14]  uses  a  labeled  Voronoi  diagreun  to  cheurac- 
terize  geometric  relations  among  line  segments.  As  will  be 
described  in  Section  2,  this  labeled  Voronoi  diagram  ceun  be  used 
for  various  picture  processing  tasks. 

The  algorithm  given  in  [14]  is  based  on  the  continuous 
Euclidean  plane.  The  objective  of  this  paper  is  to  describe  aui 
algorithm  to  compute  a  labeled  Voronoi  diagreun  in  the  digital 
plane  where  a  set  of  digital  straight  line  segments  is  given  as 
input. 

Our  algorithm  uses  a  specialized  dist£Ufice  transform  in  the 
digital  pl2me,  as  opposed  to  the  usual  distance  transforms  in 
digital  image  processing  >  i.e.,  city-block  [15] ,  chessboard,  or 
octagonal  (a  combination  of  city-block  and  chessboard,  which  is 
generated  by  using  them  alternately  at  even  euid  odd  time  steps) 
distance  transforms. 

The  main  problem  with  these  distance  transforms  is  their 
lack  of  invariance  under  rotation.  This  is  illustrated  in  Figure 
1.  Figure  la  illustrates  the  effect  of  rotation  on  the  voronoi 
edges  conqputed  using  the  city  block  distemce;  Figure  lb  is  a  simi¬ 
lar  illustration  using  the  chessboard  distance.  The  city-block, 
chessboard,  and  octagonal  distance  tremsforms,  as  well  as  the  other 
digital  distance  transforms,  in  general  also  fail  to  produce  con¬ 
nected  medial  euces. 


A  Euclide2ui  digital  distance  treuisform  is  needed  in  computing 
point  distances  in  the  digital  plctne  in  order  to  solve  the  problem 
of  the  invarieuice  under  rotation  euid  the  problem  of  the  connec¬ 
tivity  of  medial  cuces.  Monteuiari  [16]  proposed  a  quasi-Euclideeui 
distcuice  transform  for  obtaining  a  connected  skeleton  on  the  digi¬ 
tal  space.  Recently,  Danielsson  [17J  gave  an  algorithm  for  com¬ 
puting  a  digital  Euclidean  transform  of  a  given  binary  image. 

Figure  2  shows  the  Voronoi  boundary  resulting  from  this  digital 
Euclidean  transform,  for  the  same  input  as  used  in  Figure  1.  As 
we  Ceui  see,  the  shape  of  the  resulting  voronoi  boundary  is  in¬ 
variant  under  rotation. 

We  shall  define  a  special  labeling  scheme  for  the  parts  of 
line  segments.  We  shall  then  show  that  if  we  use  a  suitably 
chosen  distcuice  transfom  to  define  the  Voronoi  diagram  of  a  set 
of  labeled  digital  line  segments,  we  will  be  able  to  partition 
the  pixels  on  the  edges  of  the  Voronoi  polygons  based  on  the 
l£ibels  of  the  neighbors  of  the  voronoi  edges.  The  leUseling  scheme 
and  partitioning  algorithm  are  discussed  in  Section  2. 

In  Section  3,  a  general  algorithm  for  computing  the  Icd^eled 
Voronoi  diagreun  of  a  set  of  l6d}eled  lines  or  curves  in  a  given 
digital  image  is  presented.  In  the  algorithm,  we  first  construct 
the  Voronoi  diagram  for  the  labeled  segments;  then  each  pixel  which 
is  on  the  edge  of  a  Voronoi  polygon  obtains  a  special  leUsel  depend¬ 
ing  on  the  labels  of  the  pixels  in  its  neiahborhood. 


.  Ledaeling  and  Partitionin 


The  structure  of  Voronoi  dxagrauns  has  found  wide  applica¬ 
tion  in  euieU-yzing  the  shapes  or  the  patterns  of  a  given  set 
of  objects  in  either  the  Euclidean  or  digital  plane.  The  area 
of  a  Voronoi  region  C2ui  be  used  to  characterize  the  distance 
of  an  object  from  its  neighbors,  euid  the  shape  of  the  region 
C2U1  be  used  to  determine  the  number  of  neighbors. 

Lee  [5]  has  used  the  structure  of  the  Voronoi  diagram  as 
a  method  of  confuting  the  medial  axis  of  a  given  polygon.  The 
bisector  defined  by  two  objects  is  the  locus  of  all  points  having 
equal  disteuice  from  the  two  objects.  Thus  for  a  set  of  two 
objects  the  bisector  £uid  the  Voronoi  edge  are  identical,  auid 
for  a  set  of  three  or  more,  all  Voronoi  edges  are  segments  of  bi¬ 
sectors.  Lee  shows  that  the  gecxnetric  properties  of  bisectors 
may  be  used  for  extracting  the  medial  axis  of  a  polygon  from  its 
Voronoi  diagreun  .  In  fact,  the  medial  axis  is  a  subset  of  the 
Voronoi  diagreua  in  the  case  of  a  simple  polygon.  Matsuyama  [14] 
also  studied  the  geometric  properties  of  Voronoi  edges  and  sug¬ 
gested  a  method  for  partitioning  and  labeling  Voronoi  edges. 

Suppose  we  are  given  a  pair  of  straight  line  segments  in  the 
Euclidean  plane.  In  general,  their  Voronoi  diagreun  is  as  in 
Figtire  3,  where  the  space  is  divided  into  two  half  pl£Uies 
(Voronoi  regions)  by  a  continuous  curve  (bisector  or  Voronoi 
edge) .  The  Voronoi  edge  is  con^osed  of  five  segments  at  most. 
Using  notation  similar  to  that  in  [18],  these  segments  are 
fi(a,c),  B((a,b),c),  B ( (a,b) , (c,d) ) ,  B((a,b),d),  and  B(b,d). 


Here  a,  b,  c,  and  d  denote  the  endpoints  of  the  line  segments, 
and  (a,b) ,  (c,d)  are  their  interiors,  which  are  open  line  seg¬ 
ments.  B(p,g)  denotes  the  bisector  of  the  two  elements  p  and 
g.  B (endpoint, endpoint)  is  a  straight  line  perpendicular  to 
the  straight  line  containing  these  two  endpoints,  B (open-line- 
segment,  open-line-segment)  is  the  bisector  of  the  angle  defined 
by  the  two  line  segments,  and  B (endpoint, open-line-segment)  is 
a  parabola  whose  focus  is  the  endpoint.  By  6(Prq)  we  denote 
the  part  of  the  Voronoi  edge  arising  from  elements  p,q.  Lee 
[5]  discusses  several  geometric  properties  of  the  Voronoi  dia- 
grcun  of  a  set  of  straight  line  segments. 

As  we  can  see  from  the  example  in  Figure  3,  the  Voronoi  dia- 
greun  represents  various  geometric  properties  of  the  position  and 
length  of  the  given  line  segments  which  may  be  used  for  shape 
emalysis.  But  since  each  Voronoi  edge  is  composed  of  several 
segments  representing  different  parts  of  the  given  line  segments, 
it  is  convenient  to  regard  it  not  as  an  individual  unit  for  shape 
euialysis,  but  as  a  collection  of  segments. 

In  general,  Voronoi  edges  can  be  partitioned  into  three  dif¬ 
ferent  classes  of  segments  by  using  geometric  leUsels  for  the 
segments  B  (endpoint, endpoint) ,  B  (endpoint, open-line-segment) , 
and  B  (open-line-segment, open-line-segment) .  As  illustrated  in 
Figure  3,  a  Voronoi  edge  changes  its  leU^el  when  it  crosses  over 
one  of  the  straight  lines  which  are  perpendicular  to  the  given 
line  segments,  and  which  cross  these  line  segments  at  their  end¬ 
points  . 

Suppose  we  are  given  a  pair  of  digitized  straight  line  seg¬ 
ments.  In  the  conventional  algorithm  for  confuting  the  connected 


Voronoi  edge  of  these  two  line  segments,  first  a  unique  label 
is  assigned  to  all  the  pixels  of  each  line  segment.  Then  the 
l2d>els  are  propagated  using  the  given  distance  trauisform. 

The  Voronoi  edges  consist  of  those  pixels  where  different 
led)els  meet.  (See  Figure  4.) 

Assume  we  want  to  label  each  voronoi  edge  pixel  according 
to  its  class  -  i.e.,  assign  it  to  one  of  the  three  classes  de¬ 
fined  eJsove.  This  Icdjeling  is  impossible  without  further  analy¬ 
sis  if  we  do  not  discriminate  the  endpoints  from  the  interior 
pixels  in  each  line  segment. 

In  our  approach,  first  we  label  the  endpoints  of  a  given 
line  segment  differently  from  all  other  pixels  interior  to  the 
line  segment.  The  classification  of  a  Voronoi  edge  pixel  then 
becomes  simple  and  can  be  done  at  the  stage  of  Voronoi  edge 
extraction.  The  details  of  labeling  line  segments  euid  classifying 
Voronoi  edges  will  be  given  in  the  following  subsections.  The 
general  method  for  Icd^eling  (label  propagation)  is  similar  to 
that  used  in  the  conventional  algorithm  described  above;  the 
results  depend  upon  the  chosen  distance  tremsform.  See  Figure  5. 


As  input,  a  binary  image  I  [M,NJ  of  size  MxN  is  given,  which 
contains  a  set  of  line  segments  (composed  of  pixels  with  value  1) 


on  a  background  of  pixels  with  value  0.  As  output,  an  image 
L[M,N]  of  size  MXN  has  to  be  computed  containing  the  labeled 
segments.  Our  procedure  can  be  described  as  follows: 
a)  The  procedure  for  initial  labeling: 

Step  1:  Apply  a  connected  component  algorithm  [15]  to  image 
I,  and  give  the  pixels  in  each  component  a  unique  label. 

Step  2:  For  all  end-of-line  pixels,  add  a  large  constant 
to  their  labels.  The  constant  must  be  large  enough  that  labels 
from  the  result  of  adding  this  constamt  will  not  be  confused 
with  other  labels  in  the  image. 

This  procedure  can  be  modified  to  handle  the  case  where  the 
input  image  I  contains  digital  curves.  In  this  case,  sharp 
corners  will  be  treated  as  end-of-line  pixels. 


2.2  Classification  of  Voronoi  edae  points 


Suppose  we  are  given  a  set  of  digital  straignt  line  segments 
which  are  labeled  according  to  the  scheme  in  Section  2.1.  Now, 
a  certain  disteuice  transform  algorithm  (described  in  Section  3) 
is  applied  to  these  labeled  line  segments  for  label  propagation. 
We  shall  call  the  image  resulting  from  the  application  of  this 
transform  the  computed  labeled  distance  transform  image  P[M,N]. 
The  pixels  in  a  computed  labeled  distance  treuisform  image  have 
certain  labels  denoting  either  an  "open-line-segment"  pixel 
or  an  "endpoint"  pixel.  The  procedure  for  transforming  image 
L[M,N]  into  the  computed  labeled  distance  treuisform  image 
P[M,N]  will  be  considered  in  Section  3. 

For  these  labels  in  the  confuted  labeled  distcuice  transform 
image  we  use  the  following  definitions.  Let  LL  be  the  B (open- 
line-segment,  open-line-segment)  pixel  lcd>el,  PL  the  B (endpoint, 
open-line-segment)  pixel  label,  PP  the  B (endpoint, endpoint) 
pixel  led>el,  and  BB  the  label  for  a  branching  point  on  a  Voronoi 
edge.  Furthermore,  let  BG  be  the  background  pixel  label,  and 
OJ  the  pixel  label  for  the  original  line  segments. 

Now,  for  classifying  an  arbitrary  pixel  in  the  computed 
labeled  distance  treuisform  image  P,  the  following  definitions 
are  used: 

1)  Two  labels  11,12  in  a  computed  labeled  distance  trans¬ 
form  image  are  said  to  be  equivalent,  11  ^  12,  iff 
2.1  and  12  belong  to  the  seune  line  segment;  i.e.,  an 
endpoint  and  an  interior  pixel  from  the  same  line  seg¬ 
ment  are  equivalent. 


2)  A  pixel  in  a  computed  labeled  distance  transform  image 
is  a  Voronoi  edge  pixel  iff  there  are  two  or  more  non¬ 
equivalent  labels  in  its  neighborhood,  where  the  defi¬ 
nition  of  neighborhood  depends  upon  the  chosen  distauice 
transform. 

3)  A  pixel  in  a  computed  labeled  distance  transform  image 
is  a  BB  pixel  iff  there  are  more  than  two  non-equivalent 
labels  in  its  neighborhood. 

4)  A  pixel  in  a  computed  labeled  distance  transform  image  is 
an  LL  pixel  iff  there  are  exactly  two  non-equivalent 
open-line-segment  labels  in  its  neighborhood. 

5)  A  pixel  in  a  computed  labeled  distance  transform  image 
is  a  PP  pixel  iff  there  are  exactly  two  non-eqai valent 
endpoint  labels  in  its  neighborhood. 

6)  A  pixel  in  a  computed  labeled  distance  transform  image  is 
a  PL  pixel  iff  there  are  one  open-line-segment  label  8.1 
and  one  endpoint  label  12,  llfill,  in  its  neighborhood. 

7)  A  pixel  in  a  computed  labeled  distance  transform  image 
is  a  BG  pixel  iff  it  is  not  a  Voronoi  edge  pixel,  and  if 
it  is  not  on  any  of  the  original  line  segments. 

8)  A  pixel  in  a  computed  labeled  distance  transform  image 
is  an  OJ  pixel  iff  it  is  on  one  of  the  original  line 
segments. 

The  classification  procedure  is  as  follows: 

AS  input,  a  computed  labeled  distance  transform  image  P[M,N] 
is  given,  which  was  obtained  by  applying  the  given  distance 


transform  to  the  image  L[M,N].  The  image  L[M,N]  was  computed 
according  to  the  procedure  in  Section  2.1. 

As  output,  a  labeled  Voronoi  diagraun  V[M,NJ  will  be  deter¬ 
mined.  Our  procedure  can  be  described  as  follows: 

For  each  pixel  P[i,j]  in  P  do: 

Check  the  neighborhood  of  P[i,j]  and  assign  the  proper 
label  to  V[i,j]  according  to  the  definitions  of  Voronoi 
edge  classification  given  above. 

In  Figure  8,  we  show  the  labeled  Voronoi  edge  of  the  pair 
of  straight  line  segments  that  was  used  in  Figure  5.  In  the 
next  section,  we  shall  describe  various  possibilities  for 
using  distance  treinsforms  for  computing  the  image  P[M,N]  from 
the  image  L[M,N]. 


3.  The  labeled  distance  transform  image 

Informally,  the  process  of  generating  the  Voronoi  diagreua  of 
a  set  of  objects  in  the  Euclideeui  plane  cam  be  described  by  a 
"grassfire  procedure"  (like  that  used  for  computing  the  medial 
axis  [19])  which  operates  as  follows.  Imagine  the  space  is  filled 
with  some  flammable  material  cind  a  fire  is  started  at  each  ob¬ 
ject  (line  segment, point , curve, ...) .  If  eui  object  consists  of 
more  than  a  single  point,  the  fire  is  started  at  all  points  on 
the  object  at  once.  We  may  asstame  each  object  has  its  own  color 
of  fire.  The  fire  spreads  in  all  directions  with  equal  speed. 

The  locations  where  fire  fronts  from  different  objects  meet 
together  form  the  Voronoi  diagraun  of  the  objects. 

The  computation  of  a  digital  Voronoi  diagram  is  closely  re¬ 
lated  to  that  of  the  medial  axis  transform.  The  main  differences 


are  that  the  medial  axis  transform  is  usually  defined  for  a 
closed  boundary  of  an  object,  and  that  the  points  on  the  medial 


axis  are  leOseled  with  certain  distance  values  given  by  the  treuis- 
form.  The  Voronoi  diagr^un,  on  the  other  heuid,  is  defined  for  a 
set  of  disjoint  objects,  and  the  points  on  the  voronoi  edges  are 
normally  not  ledseled  with  distance  values. 

In  the  two-dimensional  digital  pleme,  the  Voronoi  diagram 
of  a  given  digital  image  may  be  conputed  as  follows: 

1)  Labeling  :  Give  each  object  a  unique  label 

2)  Label  propagation:  Propagate  all  leUaeled  pixels  on  the 

object  boundaries  in  predetermined 
directions  until  different  labels 


meet 


3)  Edge  extraction:  Extract  the  Voronoi  edges  at  the 

places  where  different  labels 
meet  together. 

Normally,  in  this  procedure  either  the  city-block  (4- 
direction)  or  the  chessboard  (8-direction)  distance  trauisform 
is  used  in  lcd>el  propagation  (step  2)  .  The  well-known  problem 
when  we  use  these  distemce  treuisforms  is  non-invariauice  under 
rotation,  as  we  saw  in  Figure  1.  in  Figures  6a,  7a,  euid  8a, 
the  Voronoi  label  expeuisions  of  a  pair  of  straight  line  seg- 
segments  are  shown  after  four  iterations  using  three  differ¬ 
ent  distamce  transforms.  The  final  expemsions  are  shown  in 
Figures  6b,  7b,  and  8b. 

As  we  c^m  see  from  these  Figures,  the  labels  expemded  from 
the  objects  form  either  box-like  or  diamond-like  shapes.  These 
shapes  are  quite  different  from  the  circular  shape  that  would 
be  obtained  by  using  Euclideem  distance.  This  problem  is  less 
significant  if  octagonal  distauice  is  used  in  the  computation  of 
ledsel  propagation,  but  the  octagonal  expansion  is  still  signi¬ 
ficantly  different  from  the  Euclideeui  expemsion. 

In  the  grassfire  procedure  as  we  have  informally  defined 
it  for  the  Voronoi  diagram,  the  fire  should  spread  from  the  ob¬ 
ject  in  all  directions  with  equal  speed.  Thus  even  in  a  digi¬ 
tal  image,  we  want  to  define  a  propagation  process  in  which 
the  pixels  on  objects  expand  their  labels  with  equal  speed,  i.e., 
expemd  one  unit  distemce  in  every  direction  at  each  step.  Next, 
we  will  describe  a  parallel  algorithm  that  expemds  labels  in  this 


3.1  The  p^trallel  algorithm 


The  labeled  digital  Euclidean  disteuice  trcuisform  algorithm 
(LED)  takes  a  set  of  l^d^eled  line  segments  as  input  image; 
these  line  segments  are  given  in  image  L.  The  line  segments 
in  L  are  led^eled  according  to  our  procedure  in  Section  2.1.  At 
the  beginning  of  the  computation/  LED  trzmslates  the  input  image 
L  into  a  vector  image/  V/  as  follows, 
a)  Initial  step: 

A  pixel  L(i,j)  in  L  will  be  represented  as  a  vector  of 
three  components  in  V,  V(i/ j)  =  (X/y/i)  /  where  X/  y,  and  i.  are 
integers.  The  Euclideeui  distance/  d=sqrt{x*x  +  y*y) /  is  the 
distance  from  pixel  L(i/j)  to  its  nearest  (according  to  the 
Euclidean  distance  d)  object/  say  L(h/k)  ,  such  that  [i-'hl  is 
equal  to  x,  [k-jj  is  equal  to  y,  and  L(h,k)  has  label  i.  Ini¬ 
tially/  V(i/j)  =*  (0/0/L(i/j))  iff  L(i/j)  is  not  zero  (not  a 
background  pixel)  and  V(i/j)  »  ($/$/0)  if  L(i/j)  is  zero,  where 
$  stands  for  a  very  large  number. 

In  the  algorithm  we  use  am  expansion  counter  W_count 
to  control  the  label  expansion  speed.  W_count  is  set  to  zero 
at  the  beginning.  Without  this  expeuision  counter,  the  distances 
of  the  l2U3el  propagations  would  not  be  correct  emd  the  resulting 
Voronoi  edges  would  have  the  seune  shapes  as  if  we  used  the  chess¬ 
board  distance  tremsform.  The  led>el  expeuision  resulting  from  not 
using  W_count  is  illustrated  in  Figures  Sa  and  9b. 

For  each  parallel  step  of  the  label  expansion  procedure, 
W_count  will  be  incremented  by  1,  and  each  V($/$/0)  will  be  modi¬ 
fied  as  follows: 


b)  LeUsel  expeuxsion  step  (while  any  V($,$,0)  remains  in  V)  : 

Let  VO  be  the  pixel  V(i,j)  with  vector  V(xO,yO,iO) ; 

let  VI,...,  V8  with  V(xl,yl,ill}  , . . .  ,V(x8,y8,i8)  be  its  3x3  neighbors 
V5  VI  V6  (x5,y5,i5)  (xl,yl,ill)  (x6,y6,il6) 

V3  VO  V4  =  (x3,y3,i3)  (xO,yO,JlO)  (x4,y4,14) 

V7  V2  V8  (x7,y7,i7)  (x2,y2,Jl2)  (x8,y8,i8) 

In  the  algorithm  we  confute  the  values,  val(x,y)a(x*x)  + 
(y*y) t  for  all  of  the  following  eight  pairs  (x,y) ; 

(x5+l,y5+l)  (xl+l,yl)  (x6+l,y6+l) 

(x3,  y3+l)  (x4,  y4+l) 

{x7+l,y7+l)  (x2+l,y2)  (x8+l,y8+l) 
and  let  k  be  such  that  val(x,y)  is  minimal  for  the  corresponding 
pair  (x,y)  in  conparison  to  all  the  other  7  computed  values.  (If 
k  is  not  uniquely  defined  in  this  way  then  we  select  k  according 
to  larger  x-component,  e.g.)  If  (x,y)  is  the  pair  giving  rise 
to  k,  and  if  SQRT(val(x,y) )  is  less  or  equal  to  W_covint  +  3QRT 
(2)-l,  then  the  center  pixel,  V(i,j),  will  obtain  the  value  (x, 
y,)lk);  otherwise  there  will  no  change  for  V(i,j)  .  For  exan^le, 
if  x  is  x5+l  and  y  is  y5+l,  then  V(i,j)  will  be  V(x5+l,y5+l,i5) . 

If  X  is  xl+1  and  y  is  yl,  then  V(i,j)  will  be  V(xl+l,yl,Jll)  . 

This  parallel  ladxel  propagation  stops  when  there  is  no  value 
V($,$,0)  remaining  in  the  image  V,  i.e.,  it  stops  after  max{M,N} 
steps  at  most. 

c)  Edge  extraction  step: 

The  Voronoi  edges  can  be  meurked  by  checking  the  leU^el 
changes  in  the  third  component  of  each  V(i, j) ;  see  Section  2.2. 


The  labeled  distance  transform  image  resulting  from  this  algo¬ 
rithm  ceui  be  used  for  meuiy  purposes.  The  Voronoi  diagram  of  input 
image  I  can  be  obtained  by  outputting  all  pixels  marked  as  the 
Voronoi  edge  pixels.  Furthermore,  the  Voronoi  diagram  now  con¬ 
tains  information  about  the  distances  to  the  given  objects  which 
may  be  useful  for  reconstruction  or  shape  euialysis.  The  main 
purpose  of  this  algorithm  is  to  m2dce  it  possible  to  classify  the 
Voronoi  edges  into  the  three  classes  described  in  Section  2. 

The  vector  representation  for  this  algorithm  in  image  V  is 
a  modification  of  the  two  con^onent  vectors  which  were  used  in  [17] 
Jsy  Oeuiielsson  for  the  confutation  of  the  "Euclidean  distamce 
mapping."  For  sequential  computation,  Danielsson  used  a  vector 
of  two  components  for  computing  distance  transforms  in  binary 
digital  images  in  two  passes,  with  two  picture  scans  in  each  pass 
(first  pass  from  top  to  bottom,  second  pass  from  bottom  to  top) . 

Our  algorithm  can  also  be  implemented  in  two  passes.  This  will  be  de 
scribed  in  Section  3.2. 

In  Figures  8a  and  8b,  an  exanfle  of  our  algorithm  is  shown  (for 
four  iterations  and  for  the  final  expansion)  .  Note  that  when  o\ir 
algorithm  is  performed  in  parallel,  it  has  time  complexity  0(max 
{M,N})  on  a  mesh-connected  array  processor. 


3.2  The  sequential  algorithm 

The  steps  in  the  sequential  algorithm  are  described  in  the 
seune  order  as  we  used  for  the  parallel  algorithm  in  Section  3.1. 
An  inportemt  feature  of  the  sequential  algorithm  in  comparison 
to  the  parallel  algorithm  is  that  we  do  not  use  a  step  counter. 

a)  Ini  tied,  step: 

The  sequential  algorithm  translates  the  input  image  L 
into  the  vector  image  V  and  initializes  each  vector  V(i,j)  just 
as  the  parallel  algorithm  does. 

b)  First  picture  pass: 

Let  VO  be  the  pixel  v(i,j)  with  vector  V(xO,yO ,i0) ; 
let  VI,..., V5  with  V(xl,yl,J.l)  , . . .  ,V(x5,y5,il5)  be  its  2x3  neigh¬ 
bors: 

V4  VI  V5  (x4,y4,i4)  (xl,yl,il)  (x5,y5,i5) 

V2  VO  V3  (x2,y2,i2)  (xO,yO,ilO)  (x3,y3,i3) 

1)  Scan  from  left  to  right: 

We  compute  the  values  val(x,y)=(x*x)+(y*y)  for  all  of 
the  following  five  pairs  (x,y) : 

(x4+l,y4+l)  (xl+l,yl)  (x5+l,y5+l) 

(x2,  y2+l)  (xO,  yO) 

Let  k  be  such  that  val(x,y)  is  minimal.  Then  pixel  V(i,j) 
gets  the  value  (x,y,ik). 

2)  Scan  from  right  to  left: 

This  is  analogous  to  the  left-to-right  scan  for  the  fol¬ 
lowing  five  pairs  (x,y) : 

(x4+l,y4+l)  (xl+l,yl)  (x5+l,y5+l) 

(xO,  yO)  (x3,  y3+l) 


emd  with  the  same  replacement  rule  for  V(i,j) 


c)  Second  picture  pass: 

Let  VO  be  pixel  V(i,j)  with  vector  V(xO,yO,ilO)  ,  and 
let  VI,...,  V5  with  V(xl,yl,£l) , . . .  ,V(x5,y5,£.5)  be  its  2x3 
neighbors  in  the  following  positions: 

V2  VO  V3  _  (x2,y2,i2)  (xO,yO,JlO)  (x3,y3,il3) 

V4  VI  V5  (x4,y4,il4)  (xl,yl,il)  (x5,y5,il5) 

1)  Scan  from  right  to  left,  using  the  pairs 

(xO,  yO)  (x3,  y3+l) 

(x4+l,y4+l)  (xl+l,yl)  (x5+l,y5+l) 

2)  Seem  from  left  to  right,  using  the  pairs 

(x2,  y2+l)  (xO,  yO) 

(x4+l,y4+l)  (xl+l,yl)  (x5+l,y5+l) 

d)  Edge  extraction  step: 

The  Voronoi  edges  can  now  be  marked  by  checking  the 
l^d}el  in  the  third  component  of  each  V(i,j). 

The  serial  algorithm  for  confuting  the  ledjeled  Voronoi 
diagreun  has  time  complexity  O(MxN) ,  altogether.  In  Figure  10, 
the  application  of  the  algorithm  is  illustrated  for  an  input 
image  containing  nine  digital  straight  line  segments. 


4.  Svanmary 

In  this  paper,  a  special  l£j3eling  scheme  for  a  set  of  digital 
line  segments  was  explained.  This  labeling  scheme  can  be  ex¬ 
tended  to  an  arbitrary  set  of  digital  curves.  The  ends  and  any 
sharp  comers  of  a  given  digital  curve  can  be  marked  analogously 
to  the  ends  of  digital  straight  lines.  For  example,  Fischler 
[20]  introduced  a  labeling  scheme  for  marking  the  "critical" 
points  for  a  set  of  closed  object  boundaries.  (His  l2d)eling 
scheme  was  used  to  obtain  skeletons  from  noisy  object  boundaries.) 
In  the  next  step,  by  propagation  of  the  labels  on  the  digital 
curves  we  can  partition  the  generalized  Voronoi  diagreun  of  the 
set  of  curves  according  to  their  geometric  properties.  The 
partitioned  Voronoi  edges  can  then  be  used  for  further  shape 
analysis. 

For  the  label  propagation  process,  we  gave  a  parallel  algo¬ 
rithm  for  realizing  a  Euclidean  "grassfire”  label  expemsion  in 
a  digital  image.  The  serial  implementation  of  this  parallel 
algorithm  was  also  explained.  The  resulting  sequential  algo¬ 
rithm  is  closely  related  to  the  "Euclideem  distance  mapping"  in 
[17].  The  errors  in  comparison  to  the  ideal  Euclidean  distance 
are  very  small  and  can  be  ignored  (for  the  exact  error  see  [17]). 
The  labeled  Euclidem  voronoi  diagrams  as  computed  in  this  paper 
will  be  used  for  characterizing  the  shapes  of  digital  curves  - 
e.g.,  for  conqpletion  or  reconstruction  of  broken  boundaries. 
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Figure  la,  cont'd. 


Fig.  lb,  cont'd. 
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Fig.  4b.  Analogous  to  Fig.  4a,  using  the  ohessboard 
distance  transform. 
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Fig.  5b.  Analogous  to  Fig.  5a,  using  the  chessboard 
distance  transform. 
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Pig.  10.  Labelled  Euclidean  Voronoi 
diagram  for  a  set  of  nine 
digital  straight  line  segments 


